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ABSTRACT 

The assembly of galaxies can be described by the distribution of their star formation as a function of cosmic 
time. Thanks to the WFC3 grism on HST it is now possible to measure this beyond the local Universe. Here we 
present the spatial distribution of Ha emission for a sample of 54 strongly star-forming galaxies at z ^ 1 in the 
3D-HST Treasury survey. By stacking the Ha emission we find that star formation occurred in approximately 
exponential distributions at z ^ 1, with median Sersic index of n = L0±0.2. The stacks are elongated with 
median axis ratios of b/a = 0.58 ± 0.09 in Ha consistent with (possibly thick) disks at random orientation 
angles. Keck spectra obtained for a subset of eight of the galaxies show clear evidence for rotation, with 
inclination corrected velocities of 90 to 330km/s. The most straightforward interpretation of our results is that 
star formation in strongly star-forming galaxies at z ^ 1 generally occurred in disks. The disks appear to be 
"scaled-up" versions of nearby spiral galaxies: they have EW(Ha)^100A out to the solar orbit and they have 
star formation surface densities above the threshold for driving galactic scale winds. 

Subject headings: galaxies; evolution — galaxies: formation — galaxies: high-redshift — galaxies: structure 
— galaxies: kinematics and dynamics — galaxies: star formation 



1. INTRODUCTION 

Galaxy formation is a complex process, involving star- 
bursts, mergers, and strong gas flows (e.g., Hopkins et al. 
2006; Brooks et al. 2009; Dekel et al. 2009). Furthermore, 
stellar migration and secular processes can change the struc- 
ture of galaxies at late times (e.g., Roskar et al. 2008; Grand, 
Kawata, & Cropper 2012). Therefore, even if we could per- 
fectly locate and age date every star in the Milky Way, we 
still couldn't say where and with what structural and kine- 
matic properties those stars formed. The only way to estab- 
lish where a galaxy's stars formed, and hence how it was as- 
sembled, is to map the star formation while those stars were 
forming. 

Obtaining Ha and stellar continuum maps of z > 1 galax- 
ies, with the ^ 1 kpc resolution necessary to put constraints 
on the spatial distribution of star formation, is challenging 
and has so far only been possible using a combination of 
adaptive optics and integral field units (IFUs) on 8-lOm class 
telescopes. These studies paint a complex picture: they find 
that star forming galaxies at z ^ 2 are a mix of "puffy" 
and often clumpy rotating disks, mergers, and more com- 
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pact dispersion-dominated objects (e.g., Genzel et al. 2008; 
Shapiro et al. 2008; Cresci et al. 2009; Forster Schreiber 
et al. 2009; Law et al. 2009; Jones et al. 2010; Mancini et al. 
2011). 

Some studies claim that there are trends in the structural 
properties of Ha as a function of redshift implying that the 
way galaxies assemble their stars varies fundamentally as a 
function of cosmic time. In particular, Epinat et al. (2009) and 
Kassin et al. (2012) suggest that galaxies become cooler and 
more rotation-dominated with time, z ^ 1 - being the epoch 
of "disk settling". This is interesting because most of the stars 
in the disks of galaxies like the Milky Way were formed in this 
epoch. 

It is now possible to obtain high spatial resolution (< 1 kpc) 
information on Ha emission at z ~ 1 with a high Strehl ratio, 
owing to the near-IR slitless spectroscopic capabilities pro- 
vided by the WFC3 camera on the Hubble Space Telescope 
(HST). In Nelson et al. (2012) we used data taken as part of 
the 3D-HST survey to build on previous ground-based studies 
by mapping the Ha and stellar continuum with high resolu- 
tion for a sample of 57 galaxies at z 1 and showed that star 
formation broadly follows the rest-frame optical light, but is 
slightly more extended. Here we stack the Ha maps of these 
galaxies to construct high signal-to-noise ratio (S/N) radial 
profiles and measure structural parameters of stellar contin- 
uum emission and star formation. These unique, high spatial 
resolution data from 3D-HST are combined with kinematics 
from the Near Infrai-ed Spectrometer (NIRSPEC) on the W. 
M. Keck telescope (McLean et al. 1998) We assume Hq = 70 
km s"'Mpc-', = 0.3, and D.^ = 0.7. 

2. 3D-HST AND SAMPLE SELECTION 

The 3D-HST survey, a 248-orbit Treasury program on the 
Hubble Space Telescope, supplies the 2D emission line maps 
needed to measure the spatial distribution of star formation. 
The WFC3 G141 grism provides spatially resolved spectra of 



2 



Nelson et al. 



3 



a: 

+ 1 



-1 







0.8<z<1.3 


,0 0.5 
og(SFR 


r 








3 <= 






^ o 


#- 


• 






• 

0'-- • ^'"^ 

• IT. 






• 

Fit: log(SFR) = (0.7logM.-6.8)±1cj 



log(M) [M J 



Fig. 1 . — Location of selected galaxies in the SFR-Mass plane. Six sources 
without IR photometry are not shown in the figure. The histogram shows the 
position of the selected galaxies relative to the rest of the galaxies by dividing 
a power law fit out of the distribution. The purple line shows the fit, the light 
purple region delineates ilcr. The selected galaxies appear to lie largely on 
the upper half of the 'main sequence' in SFR(IR+UV). 

all sources in the field. The wavelength range of the G141 
grism, l.lSfim < A < 1.65/im, covers the Ha emission line 
for galaxies in the redshift range 0.7 < z < 1.5. The survey 
also provides broad-band near-infrared imaging in the F140W 
filter, which samples the rest-frame R band for z ^ 1 . Using 
the same camera under the same conditions to map the both 
Ha and continuum emission allows us to compare their dis- 
tributions directly, a crucial advantage of this strategy. The 
maps have a spatial sampling of ~0.5kpc at z 1 (0''06 pix- 
els). The data presented in this paper are based on the ^ 70 
pointings, roughly half of the full data set, that were obtained 
prior to June 201 1 see Nelson et al. (2012). The survey, data 
reduction, and determination of derived quantities (e.g. red- 
shifts and masses) are described in Brammer et al. (2012) and 
van Dokkum et al. (201 1). 

The 3D-HST grism spectra have high spatial resolution 
and low spectral resolution (R ~ 130) meaning that emis- 
sion line structure reflects almost exclusively spatial structure 
(morphology) in contrast to data with high spectral resolu- 
tion where structure reflects velocity (rotation or dispersion). 
Emission line maps of galaxies are made by subtracting the 
continuum emission from the two-dimensional spectrum. The 
details of the procedure are described in (Nelson et al. 2012; 
Lundgrenet al. 2012; Schmidt et al. 2012 submitted) 

We selected galaxies with redshifts 0.8 < z < 1.3, total 
F140W AB magnitude < 21.9, and rest-frame EW(Ha) > 
100 A. We also limited the sample to galaxies with effectively 
no contaminating flux from the spectra of other nearby ob- 
jects. The position of the sample in the 0.8 < z < 1.3 star for- 
mation rate (SFR)-stellar mass plane is shown in Fig. 1. Star 
formation rates were derived from the IR -i- UV luminosities 
as described in Whitaker et al. (2012). The selected galaxies 
lie predominantly on the upper half of the 'star-forming main 
sequence' (as shown in the histogram of Fig 1) (e.g. Noeske 



et al. 2007; Daddi et al. 2007; Whitaker et al. 2012). These 
galaxies were selected in Ha not IR; while this selection may 
bias the sample, it may possibly also have the advantage that 
it adds less uncertainty to interpreting the spatial distribution 
of star formation due to dust extinction. 

3. THE SPATIAL DISTRIBUTION OF STAR 
FORMATION 

3.1. Stacked Ha and Stellar Continuum Emission 

We stack the high spatial resolution Ha maps from 3D- 
HST to create average Ha maps - increasing the S/N and 
providing for a reliable measurement of the structural prop- 
erties of Ha to large radii for this sample at z ^ 1 . We di- 
vide the sample in three bins based on their Ha siz43 in 
order to homogenize the sample and investigate trends with 
sizeOZ] The median stellar masses of the stacks are log(M*) = 
10.0, 10.2, 10.2. The galaxies are normalized by their F140W 
flux and centered according to their luminosity-weighted 
F140W image centers. We masked the [S ll] AA6716,6731A 
emission in the fit. At our spectral resolution, the FWHM of a 
line is 100 A, meaning that Ha A6563A and [N ll] A6583A 
are unresolved but Ha and [S ll] are separated by ~ 3 resolu- 
tion elements. We rotated the images based on their GALFIT- 
derived F140W position angles to align them along their ma- 
jor axes, summed them, and divided the resulting stack by the 
summed masks. Finally, in order to correct for the effects of 
the PSF, the stacks were deconvolved by adding the GALFIT- 
derived PSF-deconvolved model to the residuals of the fit, a 
method developed in Szomoru et al. (2010). 

Figs. 2 and 3 show the stacked, PSF-corrected Ha emis- 
sion. The S/N in the stacks is high, and the emission is 
clearly spatially extended. Furthermore, both the rest-frame 
R-band and Ha maps are clearly elongated. We quan- 
tified the structural properties of the stacks with GALFIT 
(Peng et al. 2002) using empirical PSFs. For the contin- 
uum and Ha we find weighted mean axis ratios of b/a{R) = 
0.55 ±0.07 and b/a{\{a )=0.58±0. 09 respectively, consistent 
with expectations for randomly oriented disks with thickness 
c/a=0.26,0.33. This thickness could be attributed to disks be- 
ing intrinsically thick at z ^ 1 , heterogeneous properties of the 
galaxies being stacked, the effect of bulges, a selection biased 
toward face-on galaxies, or noise. 

3.2. Radial Profiles are Nearly Exponential 

Fig. 3 shows the radial profiles derived from the stacks. Ra- 
dial profiles of the PSF-corrected stacks were created by sum- 
ming the flux in elliptical radii. A 30% correction was applied 
to account for [N ll] (see § 4). All uncertainties were derived 
by bootstrap resampling the stacks. 

The Ha emission ranges from exponential to much less 
centrally concentrated than exponential. The stellar contin- 
uum emission in all the stacks is nearly exponential. 

The radial equivalent width profiles (right panels of Fig. 3) 
are flat to within a factor of two and EW(Ha)> lOOA for 
all stacks within LSr^. This suggests that the high global 
EW(Ha )s seen in these galaxies are not driven primarily by 
central starbursts but by strong star formation at all radii . 

The derived Sersic indices (Fig. 2) of the Ha and stellar 
continuum profiles are n(Ha) = 1.2±0.6,0.9±0.1, n(Ha) < 1 

The sizes are measured using growth curves as described in Nelson et al. 
(2012) 

" The results are similar when the bins are defined according to the 
F140W sizes. 
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Fig. 2. — Stacks of the high spatial resolution, PSF-corrected maps of Ha emission (top). Stacking was done based on Ha size; the number of galaxies 
included in each size bin is listed (Wnoc*) as is the major axis effective radius (in kpc) of each stack measured by galfit (r,). The stacks have a weighted mean 
Sersic index «(Ha) ~ 1 and axis ratio of b/a(Ha )=0.58±0.09, consistent with disks at random orientation angles. Bottom panels show corresponding radial 
profiles with dark purple - small, medium purple - mid-sized, light purple - large (Fig. 3, §3.2). The stacked Ha emission always has Sersic index n < i. 



and «(F140W) = 1 .8 ± 0.4, 1 .4 ± 0. 1 , 1 . 1 ± 0.2 with weighted 
means of n(Ha) = 1 .0 ± 0.4 and «(F140W) = 1 .4 ± 0.2. As 
shown in van Dokkum et al. (2010), structural parameters 
measured from stacks are close to the mean values for the 
individual galaxies going into the stacks. The upper limit 
of n(Ha) < 1 for the largest Ha stack reflects the fact that 
the derived n(Ha) depends somewhat on the details of the fit 
(treatment of sky, fitting region) but always n(Ha) < 1 . 

We find that the stacked spatial distribution of Ha for 
galaxies in this sample is exponential or shallower All stacks 
have n < 2, implying a bulge fraction of less than 20% (van 
Dokkum et al. 1998). The radial Ha profiles (Fig. 3, left), 
seem to show a trend toward lower n(Ha) with increasing 
size. However, this trend is not statistically significant given 
the errors in the derived Sersic indices. Additionally, note that 
the median galaxy with EW(Ha)> 100 A has a high enough 
star formation surface density to drive winds out to ^ I r^, 
(Heckman 2002). 

3.3. Effects of Dust 

A major uncertainty in the interpretation of the radial Ha 
profile is the effect of differential dust extinction: if some 
parts of the galaxies are more obscured than others (see e.g. 
Wuyts et al. 2012), the derived radial profile of Ha would not 
reflect the radial profile of star formation. To assess the im- 
portance of this effect, we estimate the extinction as a function 
of radius. We determine the extinction from the radial color 
profile. 



The color profile is determined from the combination of the 
F140W stacks with ACS F814W stacks, which can be con- 
verted into rest-frame U -V color profiles. As shown in the 
left panel of Fig. 4, for most of the radial extent of the stacks 
the color is fairly constant, but it becomes redder inside a ra- 
dius of 2 kpc (see also Szomoru et al. 2012). To estimate 
roughly how color translates into missed star formation, we 
assume that the star formation that is not captured by Ha will 
be captured by the IR (Kennicutt et al. 2009). Using photome- 
try from Skelton et al. (2013) in prep, in addition to rest-frame 
U-V colors and IR-based SFRs from the NEWFIRM Medium 
Band Survey (NMBS, Whitaker et al. 2012), empirical rela- 
tions were derived for the translation of rest-frame U-V colors 
into star formation corrections: 

log(SFR(IRH-Ha))-log(SFR(Ha)) - 0.3 x (U-V) (1) 

The dust-corrected radial profiles are shown in the right 
panel of Fig 4. The smallest galaxies appear to have star 
formation which is marginally steeper than exponential. The 
large galaxies have the largest radial gradients but their im- 
plied star formation is still less steep than exponential. Al- 
though this analysis assigns the entire observed color gradient 
to extinction, the centers of the radial color profiles could be 
red because of dust or age. Importantly, the Sersic indices of 
the star formation in these stacks based on the implied dust 
correction remain close to one - « = 1.4, 1.1,0.6 for each of 
the stacks, weighted mean=1.0±0.4 So, even when account- 
ing for dust, the averaged star formation in these galaxies has 
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Fig. 3. — The stacks of Ha (left) and rest-frame R-band (middle) emission have nearly exponential (or shallower) radial profiles. The EW(Hci) profile 
(discussed in §3.2) is shown in the right panel. The bottom panels show the profiles normalized by their effective radius. The horizontal line is the SFR surface 
density criterion for driving large-scale outflows (Heckman 2002) As in Fig. 2, the dark, medium, and light colors correspond to the small, medium, and large 
stacks respectively as shown at the top. The spatial distribution of Hq emission is exponential or shallower 



a nearly exponential spatial distribution. 

4. KINEMATICS 

The flattening of the stacks and the exponential Ha profiles 
suggest that the star formation occurs in rotating disks. To 
test this, we measured kinematics for a subset of this sam- 
ple using NIRSPEC on the Keck II telescope on April 9-10, 
2012. The sample comprises eight galaxies chosen from the 
small (1), middle (5), and large (2) stack, which have sizes, 
masses, Sersic indices, star formation rates, and axis ratios 
representative of the sample as a whole. 

We used the low dispersion mode of NIRSPEC with 0.5" 
seeing and a slit width of 0.7", giving a spectral resolution 
of (T ^ 80km/s in the J band (compared to ^ 500km/s in the 
grism spectra). The slit was aligned along the major axis of 
each galaxy and observations were conducted in a series of 
four 900s exposures, dithering along the slit. The data reduc- 
tion followed standard procedures for long slit spectroscopy 
(see, e.g., van Dokkum et al. 2004). We extracted kine- 
matic information from the two dimensional spectra by fit- 
ting a gaussian to the Ha and [N ll] emission simultaneously 
at each spatial position along the slit. The median [N Il]/Ha 
is 0.3. Assuming the velocity shear in the two dimensional 
spectra is due to rotation, we take the rotational velocity to 



be the velocity difference between the geometrical center of 
the galaxy and the maximum velocity. We correct the rotation 
velocities for inclination angle using the 3D-HST axis ratios. 
We calculate one dimensional velocity dispersions by fitting a 
gaussian to the Ha emission in the central row of each spec- 
trum. We calculate the instrumental cr to be ^ 80km/s by fit- 
ting a gaussian to the sky lines. The intrinsic velocity disper- 
sion is calculated by subtracting the instrumental dispersion 
from the measured dispersion in quadrature. 

We find that all galaxies show velocity shear in their two 
dimensional spectra, with derived rotation velocities of 90- 
235 km/s uncorrected and 110-330 km/s corrected for incli- 
nation[3 Fig. 5 shows the rotation curves for one example 
galaxy for each stack. These spectra are of the best quality 
but appear to be representative of the eight. The velocity pro- 
file of the smallest galaxy in Fig. 3 does not show evidence for 
a turnover, which means the measured maximum velocity is a 
lower limit on the rotation velocity at large radii. We conclude 
that the kinematics are consistent with the disk interpretation 
of the structural properties of the Ha emission. Although 
different classification methods make it difficult to perform 

Measured velocity dispersions are also relatively high, but are difficult 
to interpret given the low spatial resolution of the Keck spectra. 
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Fig. 4. — The left panel shows the radial color profiles for each of the stacks with black for the small stack, dark gray for the medium stack, and light gray for 
the large stack. The right panel shows the implied dust corrected radial profiles of star formation surface density with the colors as in Fig. 2 and 3. Error bar on 
the right panel denotes a typical uncertainty. Dust-con'ected radial profiles of star formation remain nearly exponential. 



a quantitative comparison, our results are qualitatively con- 
sistent with the finding that a large fraction of star-forming 
galaxies at z ^ 1 appear to be rotating, (e.g., Wisnioski et al. 
2011; Epinatet al. 2012; Swinbanket al. 2012). 

5. DISCUSSION 

The central result of this paper is that the radial distribu- 
tion of stacked Ha emission in z ^ 1 galaxies is close to ex- 
ponential out to ^ 10 kpc. Combined with the axis ratios 
of the stacks and the kinematics of a subset of the sample, 
the most straightforward interpretation is that star formation 
seems typically to occur in disks at z 1 at least for galaxies 
with high EW(Ha). The factor of ^ 10 increase in the star 
formation rates of galaxies from z = to z = 1 (e.g., Damen 
et al. 2009; Fumagalli et al. 2012) is apparently driven by 
increased star formation activity in disks rather than a much 
greater prevalence of merger-driven central starbursts - c on- 
sistent with other studies (e.g., Rodighiero et al. 201 1; Wuyts 
etal. 2011,2012). 

This result raises a number of questions. First, it is not clear 
how these galaxies are related to z ^ galaxies. An average 
galaxy in this sample has a distribution of stars with a Ser- 
sic index of n = 1 .4. If this average galaxy forms stars with 
n = 1 .0, it will have a lower Sersic index at later times. Either 
we are witnessing the build up of only the latest of late-type 
galaxies; these galaxies substantially change where their star 
formation is occurring between z ^ 1 and z ^ 0; or their stars, 
once formed, must migrate into a different configuration. In 
other words, if these are the ancestors of typical z ^ spiral 
galaxies, their bulges need to be built in some way besides the 
star formation we are seeing. This could be accomplished by 
mergers or secular evolution (e.g., Boumaud, Elmegreen, & 
Martig 2009; Ellison et al. 2011). These observations could 
be understood in the context of the following picture: at z > 1, 
disks are gas-rich and turbulent. As redshift decreases, the 
gas fraction decreases, and a bar instability grows, increasing 
the central concentration and central vertical velocity disper- 



sion, (e.g., van den Bosch 2001; DeBuhr, Ma, & White 2012; 
Forbes, Krumholz, & Burkert 2012) 

We also note the somewhat surprising presence of a popula- 
tion of large, rapidly-rotating disks with relatively low stellar 
masses at z ^ 1. These galaxies have mean rf,(HQ;)=7.3kpc, 
v=240km/s, a = 89km/s, and = 2.2 x IO'^Mq, mean- 
ing they are disks 8 Gyr ago, larger in size than the Milky 
Way (Drimmel & Spergel 2001) and thicker (Freeman 1987), 
with ~ 1/3 of the stellar mass (McMillan 2011) and similar 
or higher maximum rotation velocities (Bovy, Hogg, & Rix 
2009). Both what these galaxies become in the local uni- 
verse and how such extended star forming disks were made 
in an epoch of high disk turbulence are open questions. The 
small disks on the other end of the size distribution are also of 
interest. With median re(HQ;)=1.9kpc, r(.(R)=2.1 kpc, = 
8.9 X IO^Mq, and stacked n(lia) = 1.2, they more closely re- 
semble compact versions of disky star formers than merger- 
driven star formation in spheroids (see also van der Wei et al. 
(2011).) 

There are several caveats. First, our sample is not com- 
plete in mass or in SFR(UV-i-IR). In future papers (using the 
full 3D-HST survey) we will study the distribution of star for- 
mation as a function of these parameters. Second, the distri- 
bution of star formation in the stacks is not necessarily repre- 
sentative of individual galaxies. In individual galaxies the star 
formation is clumpy (e.g. Genzel et al. 2008, 2011; Forster 
Schreiber et al. 2009; Nelson et al. 2012, Fig. 1) and it is 
difficult to quantify the structure. Stacking these clumpy ob- 
jects could be hazardous, but it may also provide more insight 
than can be gleaned from individual galaxies: assuming that 
the clumps are transient and "light up" a part of the underly- 
ing gas disk for a short time (as in e.g. Wuyts et al. 2012), 
our stacking technique effectively produces a time-averaged 
map of the star formation in the galaxies. Finally, dust atten- 
uation remains a key uncertainty, both in the selection and in 
the interpretation of the profiles. This can be addressed with 
maps of the IR emission at the full ALMA resolution, or by 
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Fig. 5. — Shown here are the rotation curves (middle row) derived from the NIRSPEC spectra (insets in comers), one example from each stack with the implied 
rotation velocities (vcorr - coiTected and vuncorr - uncoiTected for inclination) and velocity dispersion (cr) listed in km/s. The top row shows the coiTesponding 
rest-frame R-band images and the bottom row shows a false color image with the Ha emission in red and stellar continuum in blue. In these rows, elHpses mark 
the R-band and Ha effective radii (rj(R), reOia )) respectively and gray arrows have a scaled length corresponding to the projected specific angular momentum. 



measuring (spatially-resolved) Balmer decrements. 
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